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At low temperatures, the thermal conductivity 
of spin excitations in a magnetic insulator can ex¬ 
ceed that of phonons. However, because they are 
charge neutral, the spin waves are not expected to 
display a thermal Hall effect in a magnetic field. 
Recently, this semiclassical notion has been up¬ 
ended in quantum magnets in which the spin tex¬ 
ture has a finite chirality [Big. In the Kagome 
lattice, the chiral term generates a Berry curva¬ 
ture. This results in a thermal Hall conductivity 
K xy that is topological in origin. Here we report 
observation of a large n xy in the Kagome mag¬ 
net Cu(l-3, bdc) which orders magnetically at 1.8 
K [6-83. The observed K xy undergoes a remark¬ 
able sign-reversal with changes in temperature or 
magnetic field, associated with sign alternation 
of the Chern flux between magnon bands. We 
show that thermal Hall experiments probe inci¬ 
sively the effect of Berry curvature on heat trans¬ 
port ms®. 

In magnets with strong spin-orbit interaction, compe¬ 
tition between the Dzyaloshinskii-Moriya (DM) exchange 
D and the Heisenberg exchange J can engender canted 
spin textures with long-range order (LRO). Katsura, Na- 
gaosa and Lee (KNL) [1] predicted that, in the Kagome 
and pyrochlore lattices, the competition can lead to a 
state with extensive chirality \ = • Sj x Sj, (S; is 

the spin at site i) and a large thermal Hall effect. Subse¬ 
quently, Matsumoto and Murakami (MM) [3l 0] amended 
KNL’s calculation using the gravitational-potential ap¬ 
proach mm to relate n xy directly to the Berry cur¬ 
vature. In the boson representation of the spin Hamil¬ 
tonian, x induces a complex “hopping” integral t = 
y/J 2 + D' 2 -e xcl> with tan cf) = D/J (Fig. [lj\, inset) [HOSES]. 
Hence as they hop between sites, the bosons accumu¬ 
late the phase <f>, which implies the existence of a vector 
potential A(k) permeating k-space. The corresponding 
Berry curvature H(k) = Vk x A(k) acts like an effective 
magnetic field that imparts a transverse velocity to a 
magnon wave packet driven by an external gradient or 
force. Analogous to the intrinsic anomalous Hall effect in 
metals El, the transverse velocity appears as a thermal 
Hall conductivity n xy that is dissipationless and topolog¬ 
ical in origin, but now observable for neutral currents. A 
weak K xy was first observed by Onose et al. [2] in the py¬ 
rochlore ferromagnet LU 2 V 2 O 7 . A large K xy was recently 
reported in the frustrated pyrochlore Tb 2 Ti 207 E5]. 

Each magnon band n contributes a term to n xy with a 
sign determined by the integral of f2(k) over the Brillouin 
zone (the Chern number). Recently, Lee, Han and Lee 


(LHL) m calculated how K xy undergoes sign changes as 
the occupancy of the bands changes with T or B. Thus, 
in addition to providing a powerful way to probe the 
Berry curvature, K xy also yields valuable information on 
the magnon band occupancy. Sign changes from next- 
nearest neighbor exchange were studied in Ref. [10’. 

The Kagome magnet Cu(l,3-benzenedicarboxylate) [or 
Cu(l,3-bdc)] is comprised of stacked Kagome planes sep¬ 
arated by d = 7.97 A j6H8] (Fig. pjA.). The spin-1 Cu 2+ 
moments interact via an in-plane ferromagnetic exchange 
J = 0.6 meV. Spin-orbit coupling and the absence of in¬ 
version symmetry result in a DM exchange D = 0.09 
meV ( D/J = 0.15), with the DM vector D nearly || c [8]. 
The weak antiferromagnetic exchange between planes 
J c ~ 1 peV [8] , mediated by a benzene linker via a 0-5C- 
O bond, leads to a critical temperature Tc ~ 1.8 K. In 
the ordered state, the magnetization M lies in-plane but 
alternates in sign between planes. A weak B ~ 0.05 T 
is sufficient to align the moments [8]. Crystals (typically 
4x4x0.4mm 3 ) grow as transparent, bright blue, hexago¬ 
nal platelets with the largest face parallel to the Kagome 
planes (spanned by vectors a and b). The thermal con¬ 
ductivity tensor Ky is defined by J® = Kij(—djT), with 
the thermal current density. We apply the thermal 
gradient — VT || x || (a, b), with B || z || c. 

As we cool the sample in zero B, the thermal conduc¬ 
tivity k (nearly entirely from phonons) initially rises to 
a very broad peak at 45 K (Fig. m Below the peak, 
k decreases rapidly as the phonons freeze out. Starting 
near 10 K, the spin contribution k s becomes apparent. 
As shown in Fig. |Tp, this leads to a minimum in n near 
Tc followed by a large peak at ~ \Tc- Factoring out the 
entropy, we find that k/T (red curve) increases rapidly 
below Tc . This reflects the increased stiffening of the 
magnon bands as LRO is established. Below 800 mK, 
the increase in k/T slows to approach saturation. The 
open black circles represent the phonon conductivity K p h 
deduced from the large -B values of k xx (T, H) (see be¬ 
low). Likewise, K p h/T is plotted as open red circles. The 
difference k — n p h is the estimated thermal conductivity 
of magnons k s in zero B. 

Given that Cu(l,3-bdc) is a transparent insulator, it 
exhibits a surprisingly large thermal Hall conductivity 
(Fig. |). Above Tc, the field profile of K xy is non¬ 
monotonic, showing a positive peak at low B , followed by 
a zero-crossing at higher B (see curve at 2.78 K in Fig. 
[2j3). We refer to a positive K xy as “p-type”. Below Tc, 
an interesting change of sign is observed (curves at 1.74 
and 0.82 K). The weak hysteresis, implying a coercitive 
field <1500 Oe at the lowest temperatures, is discussed in 
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the supplementary information (SI). This sign-change is 
investigated in greater detail in Sample 3 (we plot K xy /T 
in Fig. (2tl). The curves of K, xy /T above Tq are similar to 
those iirSample 2. As we cool towards Tc, the peak field 
H p decreases rapidly, but remains resolvable below Tc 
down to 1 K (Fig. §))■ However, below 1 K, the p-type 
curve is swamped by a rising n-type contribution that 
eventually dominates below 0.6 K. Some features are in 
qualitative agreement with the calculations of LHL [§], 
but there are significant differences as well. 

To relate the thermal Hall results to magnons, we next 
examine the effect of B on the longitudinal thermal con¬ 
ductivity k xx . As shown in Fig. UK, k xx is initially 
^-independent for T > 10 K, suggesting negligible in¬ 
teraction between phonons and the spins. The increas¬ 
ingly strong B dependence observed below 4 K is high¬ 
lighted in Fig. [3j3. Despite the complicated evolution of 
the profiles, all the curves share the feature that the 13- 
dependent part is exponentially suppressed at large B, 
leaving a ^-independent “floor” which we identify with 
K ph (T) (plotted as open symbols in Fig. [lj3) . Subtract¬ 
ing the floor allows the thermal conductivity due to spins 
to be defined as n xx (T,H) = K xy (T,H) — n p h{T). The 
exponential suppression in all curves is more apparent if 
we replot the spin part as k xx /T vs. B/T (Fig. 3p). As 
may be seen, the aymptotic form at large B in all curves 
depends only on B/T. 

In the interval 0.9 K— Tc, k xx displays a H-shaped 
minimum at 13 = 0 followed by a peak at the field H p (T). 
Since k s (at 13 = 0) is falling rapidly within this in¬ 
terval due to softening of the magnon bands (see Fig. 

we associate the H-shaped profile with stiffening of 
the magnon bands by the applied B. At low enough T 
(<0.8 K), this stiffening is unimportant and the curves 
are strictly monotonic. We find that they follow the same 
universal form. To show this, we multiply each curve by a 
T-dependent scale factor s(T) and plot them on semilog 
scale in Fig. HP- In the limit of large- B, the universal 
curve follows the activated form 

K s xx ^Te-P A , (1) 

with the Zeeman gap A = gpsB where /3 = 1 /ksT, gs 


is the Bohr magneton, and g the g-factor. The inferred 
value of g (~1.6) is consistent with the Zeeman gap mea¬ 
sured in a recent neutron scattering experiment. 

For comparison, we have also plotted —n xy /T (at 0.47 
K) in Fig. HP- Within the uncertainty, it also decreases 
exponentially at large B with a slope close to A. Hence 
the exponential suppression of the magnon population 
resulting from A is evident in both k xx and n xy . This 
provides strong evidence that the observed thermal Hall 
signal is intrinsic to the magnons. 

According to MM m,«*v (in 2D) is given by 

01.2 rp _ 

Kxy = C 2(Pn)^n(k) (2) 

n,k 

where c 2 (x) = (1 + x)(logL±p) 1 2 3 4 5 - (logx) 2 - 2Li 2 (—x), 
with Li 2 (x) the polylogarithmic function, and p n = 
l/[eP( En ~^ — 1] is the Bose-Einstein distribution, with 
E n the magnon energy in band n and g the chemical 
potential. 

LHL have calculated k xv (T,B) applying the 
Holstein-Primakoff (HP) representation below and above 
Tc, and Schwinger bosons (SB) above Tc- In the ordered 
phase, the HP curves capture the sign changes observed 
in K xy (T, H ): a purely n-type curve at the lowest T and, 
closer to Tc, a sign-change induced by a p-type term. 
Moreover, the calculated curves at each T exhibit the 
high-field suppression, in agreement with Fig. HP- For 
Sample 3, the peak values of k x e are in agreement with 
the HP curves (0.04 K at T = 0.4 K; 0.2 K at 4.4 K). 
In the paramagnetic region, however, our field profiles 
disagree with the SB curves. Above Tc, n xy is observed 
to be p-type at all B whereas the SB curves are largely 
n-type apart from a small window at low B. The compar¬ 
ison suggests that the HP approach is a better predictor 
than the SB representation even above Tc- 

In addition to confirming the existence of a large K xy 
in the Kagorne magnet (and hence a finite, extensive y), 
the measured n xy can be meaningfully compared with 
detailed calculations based on the Berry curvature. For 
chiral magnets, K xy is capable of probing incisively the 
effect of the Berry curvature on transport currents. 
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FIG. 1: The in-plane thermal conductivity k (in zero B) measured in the Kagome magnet Cu(l,3-benzenedicarboxylate) 
[Cu(l,3-bdc)]. At 40-50 K, k displays a broad peak followed by a steep decrease reflecting the freezing out of phonons (Panel 
A). The spin excitation contribution becomes apparent below 2 K. The inset is a schematic of the Kagome lattice with the LRO 
chiral state pQ. The arrows on the bonds indicate the direction of advancing phase cf> = tan -1 D/J. Panel B plots k (black 
symbols) and k/T (red) for T < 4.5 K. Below the ordering temperature Tc = 1.8 K, the magnon contribution to k appears 
as a prominent peak that is very B dependent. Values of k and k/T at large B (identified with phonons) are shown as open 
symbols. 
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FIG. 2: The thermal Hall conductivity K xy measured in Cu(l,3-bdc). In Panel A, the sketch shows the spin-| Cu ions in 
adjacent Kagome planes linked by a benzene molecule in which 2 of the C-H bonds have been replaced by the carboxylate ion 
C-C (from Ref. [7]). The in-plane exchange J is mediated by the O-C-O bond while the weak interplane exchange J c goes 
through the 0-5C-0 bonds. Panel B plots the strongly non-monotonic profiles of K xy vs. B in Sample 2. The dispersive profile 
changes sign below ~ 1.7 K. The right scale gives k 2D / (k%/h) (per plane) obtained by multiplying K xy by dh/k% = 443.2 (SI 
units). Panel C and D show corresponding curves in Sample 3 (now plotted as K xy /T). Above Tc (Panel C), K xy /T is p type. 
The behavior below 1.90 K is shown in Panel D. At 1.09 K, the n-type contribution appears in weak B, and eventually changes 
K X y/T to n-type at all B. Right scale in D reports k x ®/(T kg/h). 
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FIG. 3: The effect of field B on k xx and scaling behavior at low T, for sample 3. The curves in Panel A show that the 
B-dependence of k, xx is resolved (in the range |B| < 14 T) only at T <~6.5 K. The expanded scale in Panel B shows that, 
near Tc (1.8 K), k xx has a non-monotonic profile with a V-shaped minimum at B = 0 (identified with stiffening of the magnon 
bands by the field). Below 1 K, however, k xx has a strictly monotonic profile that terminates in a sharp cusp peak as B —>0. 
At each T < Tc, the constant “floor” profile at large B is identified with K p h- Panel C shows that the complicated pattern in 
Panel B simplifies when plotted as k xx /T vs. B/T. Further, multiplying by a scaling factor s(T) collapses all the curves below 
1 K to a “universal” curve, as shown on log scale in Panel D. The slope at large B gives a Zeeman gap with g = 1.6. The Hall 
curve — K xy /T has a similar slope at large B. 








































